Using np.diagflat only along second dimension